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Abstract 

A simple integral representation is derived for the quasiclassical Green function 
of the Dirac equation in an arbitrary spherically-symmetric decreasing external field. 
The consideration is based on the use of the quasiclassical radial wave functions with 
the main contribution of large angular momenta taken into account. The Green 
function obtained is applied to the calculation of the Delbriick scattering amplitudes 
in a screened Goulomb field. 
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1 Introduction 


The most convenient way to take into account the external electromagnetic field in quan¬ 
tum electrodynamic processes is the use of the Furry representation. So, it is necessary 
to know the Green function G{f,f'\e) of the Dirac equation in this field. Unfortunately, 
the explicit forms of the Green functions are known only for the few potentials and nu¬ 
merical calculations should be exploited. For many high-energy QED processes the main 
contribution to the amplitudes is provided by large angular momenta. Therefore, one can 
use the quasiclassical approximation. In the present paper, the explicit expression of the 
quasiclassical Green function of the Dirac equation in an arbitrary spherically-symmetric 
decreasing external field is found. Previously the quasiclassical Green function of the Dirac 
equation has been obtained in |l|, for the case of the Goulomb field by summing the 
integral representation of the exact Green function 0 over large angular momenta. As 
it will be shown, to obtain the quasiclassical Green function, it is not necessary to know 
the exact one. It is sufficient to use the quasiclassical radial wave functions at large an¬ 
gular momenta. This method has been applied earlier in to derive Sommerfeld-Maue 
type wave functions P| used at the consideration of high-energy bremsstrahlung and pair 
production in a screened Goulomb field. The integral representation of the Green function 
obtained in our paper is convenient in analytic calculations of the amplitudes of different 
high-energy QED processes in the external field. To confirm this statement we calculate 
the Delbriick scattering amplitude (the elastic scattering of a photon in the external 
field via virtual electron-positron pairs) in a screened Goulomb field. 


Delbriick scattering is one of the few nonlinear QED processes which can be tested 
by experiment with high accuracy (see recent review [|^). At the present time Delbriick 
amplitudes are studied in detail in the Goulomb field exactly in the parameter Za a.t high 
photon energy a; S> m only; m is the electron mass, Z\e\ is the charge of the nucleus, 
Za = = 1/137 is the fine-structure constant, e is the electron charge, h = c = 1. The 

approaches used essentially depended on the momentum transfer A = \k 2 — ki\ { ki and 
^2 being the momenta of the incoming and outgoing photons, respectively). At A a; 
the amplitudes have been found in iBng by summing in a definite approximation the 
Feynman diagrams with an arbitrary number of photons exchanged with a Goulomb centre, 
and also in ii with the help of the quasiclassical Green function in a Goulomb field. At 
m -C A ~ a; the amplitudes of the process have been obtained in |]^, |T^, |T^ using the 
exact electron Green function in a Goulomb field P] in the limit m = 0. Many authors 
have performed the calculations for an arbitrary photon energy uj but only in the lowest- 
order Born approximation with respect to the parameter Za (the results obtained in this 
approximation are surveyed in |^). It turned out that Goulomb corrections at Za ~ 1 
and O' 3> m drastically change the result as compared to the Born approximation. 


The effect of screening is important only in the case of small momentum transfer A ~ 
1/rc -C m, where Tc is the screening radius of the atom. It is this range of momentum 
transfer that we consider in our paper. 
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2 Green function 


Let us consider the Green function of the Dirac equation in the external spherically- 
symmetric held V(r) : 


G(f,r'je) 


7 °(£ — V(r)) — 'yp — m + iO 


6{f- r '), 


( 1 ) 


where 7 ^ are Dirac matrices, p = —iV . We are interested in the calculation of the Green 
function at |£| S> m. Let us represent the function G in the form 


G{r, r' \ e) 


7 °(£ — V{r)) 


7 P + m D{r, r'\ e ), 


( 2 ) 


where the function -D(r, f' \ e) is 

"'I = (, - v{rW - 'P -W.r(r)| - + M 

Here a = 7^7 . As it is known (see |^), at high energies £ 3> m one can neglect V‘^{r) in (|^) 
and take into account only the hrst term of the expansion with respect to the commutator 
[dp, D(r)]. Making the cited expansion and using the representation 

[dp,V{r)] = ^[dp,H] , H = p^ + 2eV{r), (4) 

we get the following representation for the function D : 

, (5) 

where 

D^^\r,r'\e) = ( 6 ) 

— rn? . Thus, the problem reduces to the calculation of the quasiclassical Green 
function of the Schrodinger equation with the hamiltonian H. 

Let us introduce the impact parameter p = Ir'xrr'|. In high-energy processes the 
characteristic distances are \r — f'\ ~ K/rrd S> 1/m and p > 1/m . So, the corresponding 
angular momentum / ~ Kp 3> 1, and one can use the quasiclassical approximation. Besides, 
we consider the case p <C \r — r'\ . Hence, the angle either between f^and —r' or between 
r and r' is small. 


D{r, r'\ e) = 


1 - -(«,v+v': 


Gonsider the set of eigenfunctions of the hamiltonian H and use in (P) its complete¬ 
ness, replacing 5-function by the sum of products of eigenfunctions. Obviously, the main 
contribution to is provided by the functions of the continuous spectrum with large 
angular momentum values. The eigenfunction containing the plane wave with the 

momentum q and the outgoing spherical wave in its asymptotic form, can be represented 
as: 

-1 00 

'4’qi^ = + l)ui{r)Pi{cos'd). (7) 
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Here Pi{x) are the Legendre polynomials, 'd is the angle between vectors q and r. The set of 
eigenfunctions with the ingoing spherical waves in asymptotics leads to the same result for 
the Green function. In the quasiclassical approximation the functions ui{r) and 6i = S{l/q) 
are equal to (see Q]) : 


Ui{r) = sin (^qr — Itt/2 + f /2qr + \5{l/q) + A<h(r)^ , 

OO OO 

4>(r) = J V{Od( , i(p) = - j V (v'd+T) dC , A = £/,. 

r 0 

Taking into account the completeness relation, we get: 

dq 


( 8 ) 


f'\e) = 


K? — q^ + iO (27r)3 


( 9 ) 


Substituting into (H) and taking the integral over the angles of vector q with the help 
of well-known relation for the Legendre polynomials 


Pi{nin2)Pu{nin‘i) dVti = 


dvr 

2 / + 1 


Pi{n2n3)6ii> , 


where n* are unit vectors, we obtain: 

1 T 

f'\e) = 


dq 


2tt'^ 






T + *0 ^ 


— ^(2i+ l)ti,(r)ui(r')F,(nn') , 


( 10 ) 


n = f/r, n' = r'/r' . Using (H), we represent the product ui{r)ui{r') as follows : 

ui{r)ui{r') = ^ cos q{r — r') -|- f‘{r' — r)/2qrr' + A(<h(r) — $(r')) 

— -(—1)^ cos q{r + r') + P{r -|- r')/2qrr' + 2X5{l/q) + A (<h(r) -|- $(r')) 


( 11 ) 


If the angle 9 between the vectors n and —n' is small, then one can replace Pi{nn') by 
(—1)^Jo(/6'), where Ju{x) is the Bessel function. After this substitution it is clear that in 
sum over I the main contribution comes from the second term in dUD- For this term the 
summation with respect to I can be replaced by an integration. Let us make in (pTl) the 
exponential parametrization of the energy denominator : 


fX2 _ q2 iO 


= —i J exp[is(K^ — q‘^)]ds 
0 


Then the integrals over q and s can be taken by means of the stationary phase method. 
After simple calculations for the case under consideration 6^ -C 1, one obtains: 


f'\e) = 




Atikw' 


dlljQ{W) exp) i 


P{r + r') 


2Krr' 


f2A5(//fi)) + A(<h(r) + <F(r')) 


( 12 ) 


In this formula X = e/n. In the relativistic case A = -|-1 at e > 0 and A = —1 at £ < 0. 
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If the angle 6i = tt — 6 between vectors n and n' is small then one can replace Pi{nn') 
by Jo{Wi). In this case the main contribntion in the snm over I comes from the hrst term 
in ([TT|) . Since it doesn’t contain 6{l), it is possible to take the integral over I after the 
transformations similar to those performed at the derivation of (|^. At tt —0 -C 1 one has: 


f'\e) = 


— exp {m|r — f'\ + iAsign(r — r')(<h(r) — <I)(r'))} 


dvrlr*— f' 

Snbstitnting (^) into (|g) and (|), we hnd for the fnnction D aX 9 1 : 


D{r,f'\e) = 


dvTKrr' J 


dll 


MW) - i^-^^^^^^^5'{l/n)MW) 


X exp|* 

Here 6'{p) = dS{p)/dp. 


f{r + r')/2Krr' + 2X6(1/k) + A($(r) + ^(r')) 


(13) 


(14) 


At TT — 6 ^ -C 1 the fnnction D is of the following form: 

D{f, f'\e) = — [1 — sign(r — r'){V{r) — V{r')){d, n + n')/4K] 

X exp {m|r — r'I + iAsign(r — r')(<I>(r) — <h(r'))}/ 47 r|r — r'I . (15) 


Snbstitnting (|^ and (|T5|) into (§), we get the hnal resnlt for the qnasiclassical Green 
fnnction of the Dirac eqnation in spherically-symmetric external held. One has at 6 ^ -C 1 


G(r, r'\e) = 


;piK,(r+r') 


le 




dll exp) i 


f(r + r')/2Krr' + 2X6(1 /k) + A(<I)(r) + <^(r')) 


X 


7 °e + m-( 7 , n — n')(K + /^/ 2 Krr') 


Jq(19) + i 


l‘^(r — r')(j,n + n')/2rr' (16) 

Ji(W)/(W)], 


+16'(I/K)'y^ ^1 — ( 7 , n)('y,n') — ( 7 , n + n')m /k 
and at TT — 0 -C 1 

G(+r'\e) = + (f^ + ^lR){lM)/R] 

X exp {mi? + iAsign(r — r')(<h(r) — $(r'))} , R = f — f'. (17) 

In the Conlomb held V(r) = —Za/r^ we have 

2(5(p) + $(r) + <I)(r') = Za ln(4rr7p^) , 6'(p) = —Za/p. (18) 


Using ([TSp and (p! 6 [) , we hnd that onr resnlt for the Green fnnction in the Goulomb held 
is in agreement with that obtained in [|^ 


3 Delbriick scattering 

Let us apply the formulae obtained to the calculation of the Delbriick amplitudes in a 
screened Goulomb held. In the Thomas-Fermi model the screening radius Tc ~ (ma)~^Z~^/^ 
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. The characteristic impact parameter p ~ 1/A. If i? -C 1/A -C (i? is the radius of 
the nucleus), then the screening can be neglected and the amplitude under consideration 
coincides with that in the Coulomb held. If 1/A ~ 3> 1/m, then it is necessary to 

take screening into account. For this momentum transfer the main contribution to the 
amplitude is provided by impact parameters p from 1/m to Tc- The corresponding angular 
momenta I ~ cup S> 1 and the quasiclassical approximation is valid. 


Let an initial photon with momentum ki produce at the point n a pair of virtual 
particles which is transformed at the point r 2 into a photon with momentum ^ 2 . Then the 
uncertainty relation gives r ~ jfg — ri| ~ lojijv? + A^) for the lifetime of the virtual pair. 
Therefore, at uj/m? ^ the angles between ki,k 2 , r 2 and —ri are small. It is this energy 
range we consider further. According to the Feynman rules, in the Furry representation 
the Delbriick scattering amplitude reads 

M = 2ia / dri (ir 2 exp[f(fciFi — fc 2 A)] / deTrelG{f2,fi\u — e)eiG{fi,f2\ — e) , (19) 


where and 62 are the polarization vectors of initial and hnal photons, respectively, 
e = e^ 7 ^. It is necessary to subtract, from the integrand for M in (^), the value of this 
integrand at zero potential. In the following such a subtraction is assumed to be made 
and we perform it explicitly in the hnal result. The main contribution to the amplitude M 
arises at the integration over e from m to a; — m. Thus, A = +1 in the hrst Green function 
in ([I^ and A = —1 in the second one. Using the representation (|^) , it is convenient to 
rewrite eq. (p!9[) in the form 


M = ia dfidf 2 exp[i{kifi — k 2 r 2 )] / deTr 


{262*P 2 - t2h)D{r2,ri\uJ - e) 


( 20 ) 


X 


( 2 eip\ + eifci)iA(ri, r 2 | — e) + 2 iae 2 *ei J dr exp[i(fci — ^ 2 )^] J deTrD{f,f^e 


Here p 1 2 = The last term in (^) doesn’t contribute to the amplitude at high energy 

because it is independent of a; and depends only on momentum transfer A. However, the 
amplitude at a; 3> A is proportional to co (see, i.g., 0). Our further transformations are 

perform the differentiation and take the trace 


as follows. We substitute (|T^ into 
over 7 -matrices. It is convenient to direct the axis of the spherical coordinate system along 
ki + k 2 ■ In the small-angle approximation one has: dQi ^2 ~ di^ 2 ddi, 2 ^ 01^2 = d 9 i ^2 with 
(di, 2 , ki + /C 2 ) = 0 . The Bessel functions in the Green function depend on the vectors 61^2 
via the combination 9 = |di -|- ^ 2 ! only. Let us change over to the variables 9 = 9i + 92 and 
^ = ■'"idi — ^ 202 . After that it is easy to take the integral over d^. Further, to demonstrate 
the method of calculations, we consider the case of zero momentum transfer (^2 = ki = k), 
and present then the results of similar calculations for A ~ l/vc- 


3.1 Zero momentum transfer 


We set ki = /c 2 , ef = 62 and take the integral with respect to d9 with the help of the 
relation ([0, p. 732) 


■^i7ru/2 


dxxe^^^ Jy{ax)Jv{hx) = —-—Jy(a6/2c)exp 


-i{a^ + IP) 

4c 
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and also those obtained by differentiating this expression with respect to the parameters. 
Let ns make the snbstitntion of the variables in the integral representation of the Green 
fnnction: li = Kipi , I 2 = K 2 P 2 where Ki = , K 2 = ((a; — . After 

that it is convenient to pass from the variables ri and r 2 to s and x : ri = KiK 2 /[m‘^ujsx], 
r 2 = KiK 2 /[m‘^u}s{l — a:)]. As the resnlt, the integral over e becomes trivial and we get the 
following expression: 


M = 


2 iaujm? 


dx 


x{l — x) 


-I CXD 00 


1 + 


x(l — x) 


j J PiP2dpidp2 


0 0 


ds 

s 


( 21 ) 


X expA + pI)s - l/[sx(l - x)] I sin ^( 5 ( p 2 ) - <l(pi)) sp-iPl ), 


Here we have snbtracted from the integrand its valne at the field eqnal to zero. Let ns 
make the snbstitntion of the variables pi = , p 2 = pe"^/^ and deform the contonr of 

the integration with respect to s so that the integral is extended from zero to ioo. Then 
the integral over s can be taken nsing the relation (|]^, p. 739 ) 


/ cLoc 

exp[—a;(a^ + 6^)/2 — l/2x]Iu{ahx )— = 2Iy{a)Kjj{h) , a<b, 

X 

0 


where A (a:) and K^{x) are the modified Bessel fnnctions of the first and third kind respec¬ 
tively. As the result, we obtain 


M 

X 


Siaujm? r dx 

3 J x(l — a;) 

0 ^ ^ 

00 

J dr sin2(^(5(pe^/^) - (5(pe"^/^)^ Io{yi)Ko{y 2 ), 

0 


1 + 


a;(l — a;) 


p’^dp 


( 22 ) 


where = mpe"’"’”/^ [a:(l — a:)] . 

We divide the integral over r into two parts: from 0 to tq and from tq to cx) , where 
1 3> To S> 1/ {mrc) . Let us begin our calculations from the second domain. In this domain 
the main contribution is given by the impact parameters p < r^, and the field can be 
considered as the Coulomb one. Integrating over x , p, and , finally, over r, we get 


. 280 : 0 ; f , chrsin^(Zor) 

^2 = / dr -^^ = 

9m2 J sh r 

TO 

_^ 28oo;(Zo)^ - iZa) + G + ln2ro - 3/2]. (23) 

9m2 

Here '0(x) = d\nT{x)/dx , C = 0.577... is the Euler constant. 

In the first domain the difference 5(pe'^/^)—(5(pe“'^/^) is small, and it is possible to expand 
the integrand with respect to this quantity. Therefore, this domain gives the contribution 
to the amplitude in the Born approximation only. It is convenient to divide the integral 
over p into two parts: from zero to po and from po to 00, where rc 3> po S> l/(mro) . In 
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the integral from zero to po fhe field can be considered again as the Coulomb one, and the 
integrals left can be easily taken. The corresponding contribution reads: 


Mil ^ [ln(mTo/3o) + C - 11/21] . 


(24) 


In the integral from po to oo one can use the asymptotics of the Bessel functions Io{x) and 
Ko{x) at large x and extend the integration over r up to inhnity. The corresponding result 

is 


Mi2 = i 


28au 

9m2 


P 


PO 




dp. 


(25) 


The sum of (|^), (|J) and (^Sf) is equal to 


M = i 


28auj{Za) 


9mP 


7 (d6\^ 41 

ln(mpo/2) + {ZaY"^ J p ( ^ j ~ Re-?/i(l - iZa) + — 

PO ^ ^ 


(26) 


At p <C Tc the integral in (^) is equal to ln(rc/po) +A, where A is some constant. Therefore, 
the amplitude M in is independent of po , and one can put, for instance, po = 2/m. 
So, we have obtained the hnal result for the forward Delbriick scattering amplitude for an 
arbitrary screened potential. The explicit value of the constant depends on the form of the 
potential. 


Let us consider the case of the Moliere potential , which approximates the potential 
in the Thomas-Fermi model: 


Za 


r(r) =-1: 






2 = 1 


where ai = 0.1 , 0:2 = 0.55 , 0:3 = 0.35, Y = mZR^6j/121 , hi 
The corresponding scattering phase is equal to 

5(p) = ZaY^ aiKYYp). 

2 = 1 


(27) 

6 , 62 = 1.2 , 53 = 0.3 . 

(28) 


Substituting this expression into (^), we get the hnal result for the forward Delbriick 
scattering amplitude in the Moliere potential: 


M = i 


.28au:{Za) 


9m^ 


ln(183Z-R3) - C - Rep(l 


iZa) — — 
^ 42 


(29) 


As known, the imaginary part of the forward scattering amplitude of the photon is con¬ 
nected with the total cross section a of electron-positron pair production by the relation 
a = ImM/o;. Due to this relation, our formula (^) is in agreement with the result of 
for the total cross section of pair production in a screened potential. Note that the real 
part of the amplitude ( ^9]) in the screened Coulomb potential is equal to zero in contrast 
to the case of pure Coulomb potential [0, 0. 



















3.2 Non-zero momentum transfer 


At non-zero momentum transfer it is convenient to carry out the calculation in terms of 
helicity amplitudes. One can choose the polarization vectors in the form 


eA 


= ([A X z7i,2] ± i\)/V ^, A = [zTi X V2]/\[vi x V2\ \ , 


(30) 


where vi ^2 = ^ 1 , 2 !^ ■ There exist two independent amplitudes: = M and = 

M~^. In terms of linear polarization, by virtue of parity conservation, the amplitude differs 
from zero only when the polarization vectors of the initial and hnal photon both lie in the 
scattering plane (MH) or are perpendicular to it (M^). These types of amplitudes are 
related via 

Mil = M++ + M+- , = M++ - M+- . 

At A = 0 the amplitude M^~ vanishes by virtue of the conservation of the angular mo¬ 
mentum projection along the direction of motion of the initial photon, and the amplitude 
coincides with . Similar to the case of zero momentum transfer, we divide the 
integral over r into two parts: from 0 to Tq and from tq to cx) , where 1 3> Tq S> l/(mrc). 
The angle 9q between ki and ^2 is 6^0 = A/o; -C m/uj . In the domain from tq to 00 the 
held concides with the Coulomb one and the angle 9q can be neglected. The contribution 
of this domain to the amplitude coincides with M 2 (|^) , and the contribution to the 
amplitude M^~ is equal to zero. In the domain from zero to tq we split the integral over 
p into two parts again: from zero to po and from po to cx), where Tc ^ po ^ l/(mro). In 
the integral from zero to po the held can be treated as a Coulomb one and the angle 9q 
can be neglected again. The corresponding contribution to M~^~^ coincides with Mu (^) 
, and the contribution to M~^~ is equal to zero. The ehect of screening is essential in the 
last domain from po to 00 only. In this domain the main contribution to the integral over 
angles is given by 9 p/r ~ pm?/uj 3> m/uj 3> 9q. The argument of the Bessel functions 
in the expression for the Green function is 19 ~ ojp9 ~ {mpY 3> 1, and one can use the 
asymptotic expansion of the Bessel functions. It is necessary to keep two terms of the 
expansion due to the compensation. After that the integrals over 9 and the other variables 
can be easily taken, and we get the contribution of the domain under discussion to the 
amplitude M++ : 

00 / •. O 

. 280 : 0 ; r (85 


M++ = *■ 


9m^ 


/p(g^) Mp^)dp. 


PO 


(31) 


The sum of (0), (0) and is: 


M++ = 

9m^ 


(Za) 


-2 


00 /ar \2 ... 

J p [ y ] Mp^)dp-'Re?j{l-iZa) +— 

2/m ^ 


(32) 


Substitute into (^) and take the integral over p by means of formula (6.578(10)) [|T^ . 
Then, the hnal result for the amplitude M~^~^ in the case of Moliere potential reads: 

M** - - g + i (33) 




a-i 






u 


— 1 


ln(M -|- — 1) 
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where u = (A^ + jSj + f3j)/2f3i(3j . At A -C 1/rc the formula (^) coincides with (|^) . At 
m A 1/rc the formula (^31) turns into 


(34) 


.2Sauj{Za)'^ m 41 f 

in agreement with the result of iii- Similarly, for the amplitude M~^ , we obtain: 


M+“ = i 


Aauj 

Qm? 


P 


9p, 


J2(pA) dp. 


(35) 


Here the lower limit of the integral is replaced by zero since the domain from zero to po 
doesn’t contribute. For the Moliere potential one has 


M+“ = i 


2auj{Za) 


Qm? 

u{(3f + / 3 |) - 2f3.(3j 
Av? - 1 


1 ^ 


iPi - P]) In ^ 

Pi 


ln(M + yjv? — 1) 


(36) 


If A —0 then M+ (|36[) tends to zero. At m 3> A S> 1/rc the formula (|36D turns into 


M+- = i 


2auj{Za) 


9mP 


(37) 


which is in accordance with the result of iii- 


Thus, we have demonstrated on the example of Delbriick scattering in a screened 
Coulomb potential that the quasiclassical Green function obtained in the present paper 
can be used effectively at the consideration of high-energy QED processes in an arbitrary 
spherically-symmetric decreasing external held. 

We are grateful to V.M.Katkov and V.M.Strakhovenko for useful discussions. 
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